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ABSTRACT 


A  thorough  analysis  of  the  stresses,  strains  and  displacements 
surrounding  the  crack  tip  of  a  single  edge-cracked  beam  subjected  to 
Mode  Z  tensile  loading  was  performed  by  constructing  various  two- 
dimensional  finite  element  models.  Of  particular  interest  in  the 
analysis  was  1)  determination  of  the  Mode  I  stress  intensity  factor 
in  the  elastic  models  and  2)  examination  of  the  size  and  character 
of  the  plastic  zone  in  the  elastic-perfectly  plastic  models. 

The  finite  element  program  used  in  the  investigation  was  BOPACE-3D 
(The  Boeing  Plastic  Analysis  Capability  for  Engines) .  In  order  to 
model  the  singularity  effects  at  the  crack  tip,  special  /^crack  tip 
elements'1  were  utilized  Immediately  surrounding  the  crack  tip. 

Different  crack  tip  elements  were  used  for  the  elastic  models  and  the 
elastic-perfectly  plastic  models.  Both  coarse  and  fine-grid  meshes 
were  evaluated  in  the  analysis. 

The  finite  element  model  results  are  compared  to  various  other 
solutions.  Results  for  the  elastic  models  are  compared  to  available 
empirical  solutions  while  the  elastic-perfectly  plastic  results  are 
examined  with  regard  to  extent  and  character  of  the  elastic-plastic 
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CHAPTER  I 


INTRODUCTION 


1.1  General  Introduction 

Elastic  and  elastic-plastic  fracture  mechanics  problems  have 
received  considerable  interest  for  the  past  20  years.  A  great  deal 
of  work  in  elastic  fracture  mechanics  has  been  centered  around  the 
development  of  empirical  solutions  [l4,  16],  experimental 
techniques  [l7]  and  numerical  methods  [ll]  for  the  determination  of 
Mode  I  stress  intensity  factors.  Photoelasticity  is  one  of  the 
experimental  techniques  which  has  been  used  to  predict  stress  intensity 
factors  while  finite  element  modeling  is  a  numerical  method  which  made 
its  debut  in  this  area  of  fracture  mechanics  in  the  early  1970's  [llj. 

With  the  development  of  special  "crack-tip  elements,”  [2,  5,  7], 
finite  element  models  currently  represent  the  crack-tip  singularity 
effects  with  a  higher  degree  of  accuracy.  Not  only  do  the  finite 
element  models  better  represent  the  singularity  effects,  but  the  models 
are  able  to  do  so  with  coarser  meshes,  which  result  in  great  savings 
in  time  and  money. 

In  addition  to  the  development  of  elastic  crack-tip  elements, 
finite  element  programs  and  crack-tip  elements  which  represent 
plastic  behavior  have  also  been  developed  \fi\.  This  has  enabled  the 
finite  element  modeling  of  elastlc-perfectly  plastic  behavior. 
Specifically,  plastic  zone  effects  at  the  crack-tip  can  be  modeled 
for  various  loading  conditions  and  the  size  and  character  of  the 
plastic  zone  may  be  examined. 
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1.2  Purpose  of  the  Investigation 

The  purpose  of  this  Investigation  was  to  examine  and  evaluate 
recently  developed.  Improved  finite  element  techniques  which  can  be 
applied  to  fracture  mechanics  problems.  Particular  emphasis  was 
placed  on  evaluating  both  the  elastic  and  elastic-perfectly  plastic 
crack-tip  elements  in  modeling  a  single  edge-cracked  beam  subjected 
to  Mode  I  tensile  loading. 

1.3  Scope  of  the  Investigation 

This  investigation  can  be  divided  into  four  primary  parts.  The 
first  part  consists  of  generally  analyzing  the  stress,  strain  and 
displacement  distributions  surrounding  the  crack-tip  for  a  single 
edge-cracked  elastic  beam  subjected  to  Mode  I  tensile  loading.  The 
second  part  consists  of  1)  determining  elastic  Mode  I  stress  Intensity 
factors  using  the  displacements  obtained  from  the  elastic  finite 
element  model  then  2)  comparing  the  Mode  I  stress  Intensity  factors 
determined  from  the  finite  element  model  to  stress  intensity  factors 
resulting  from  a  boundary  collocation  solution  [l6].  The  third 
part  of  the  investigation  consists  of  generally  analyzing  the  stress, 
strain  and  displacement  distributions  surrounding  the  crack-tip  for 
a  single  edge-cracked  elastic-perfectly  plastic  beam  subjected  to 
Mode  I  tensile  loading.  Finally,  the  fourth  part  consists  of 
examining  the  size  and  character  of  the  plastic  zone  surrounding  the 
crack-tip  for  a  single  edge-cracked  elastic-perfectly  plastic  beam 
subjected  to  Mode  Z  tensile  loading. 


CHAPTER  II 


ELASTIC-PLASTIC  FINITE  ELEMENT  MODELING 

2.1  General  Concepts 

General  concepts  for  the  solution  of  elastic  or  elastic-plastic 
problems  by  means  of  isoparametric  finite  elements  may  be  found  in 
various  references  ^5,  6,  15 J.  Using  the  notation  found  in 
Reference  6,  the  geometry  of  an  8-node  plane  isoparametric  element 
is  mapped  into  a  normalized  square  space  (£,  n)»  through  the 
following  transformations, 

8 

x  -  £  N  (c,  n)x 
i=l  1 


and 


(1) 


8 

y  =■  z  n  U,  n)y, , 
i-1  1  1 


where 


Nt  -[(1+  eCjHl  +  nni)-(l-C2)(l  +  nn1)-(l-n2)(l  + 

+  (1-C2)(1  +  nn^d-^Jn2  /2  +  (l-n2)(l  +  «±)(l-n 2)C2  /2, 


(2) 


N^  "  the  shape  function  corresponding  to  node  1,  whose 
coordinates  are  (x.,  y.)  in  the  x-y  system  and 
(?!»  n^)  in  the  transformed  £-n  system 


Cf»  "  +1  for  comer  nodes 

0  for  raid-side  nodes. 
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The  displacements  are  Interpolated  by  the  same  shape  functions  as 
8 

u  -  E  N  (5,  n)u. 

1-1  1 

(3) 

and 

8 

v  »  E  N, (5.  n)v  . 

1-1 

The  stiffness  matrix  is  found  through  the  following  relationships. 


H  -  H 


and 


w 


3N 


_i 

3x 


0 

3N 


0 

3N 


_i 

3y 


_i 

3  y 


3N 


1 

3x  1  ’ 


(4) 


(5) 


where 


’3N1 

'3Ni' 

3x 

-  M"1 

35 

3Ni 

L  J 

3^ 

ay 

3n 

(6) 


and  the  Jacobian  matrix  £jJ  is  given  by 


’  3x  3x 

3Ni 

35  35 

"*•35 . 

3x  3^ 

■ 

3n  3n 

3n 

xi  yt 


(7) 
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The  stress  is  given  by 


H  •  H  H 


where  £dJ  is  the  stress-strain  matrix.  The  element  stiffness  matrix 
Jjcj  is  then. 


w  -  fi  fi  i»r  w  w  det | j|d£dn. 


(9) 


The  above  derivation  considers  the  elastic  case,  however,  a  similar 
derivation  at  an  incremental  level  can  be  performed  for  the  plastic 
case. 

In  modeling  fracture  mechanics  problems,  a  singularity  must 
occur  at  the  crack  tip.  Therefore,  to  obtain  a  singular  element 
which  will  be  used  at  the  crack-tip,  the  stress  and  strain  in 
equations  (8)  and  (4)  must  be  singular.  The  formation  of  crack-tip 
elements  which  model  the  crack-tip  singularity  is  the  topic  of  the 
following  section. 

2.2  Elastic  and  Elastic-Perfectly  Plastic  Crack-Tip  Elements 

One  method  used  to  create  the  singularity  in  the  quadrilateral 
Isoparametric  elements,  is  to  collapse  an  edge  of  the  element  and 
move  the  mid-side  nodes  to  the  quarter-points  as  shorn  in  Figure  1. 
For  the  element  shown  in  Figure  1,  we  have  the  following  nodal 
coordinates, 

X1  “  x7  "  *8  '  °’ 
h 

x2  "  X6  "  4 


» 


h/4'  3h/4 


and  y3  -  -y^  -  -i. 

Substituting  the  nodal  conditions  (10)  in  equations  (1)  and 
collecting  terms  yields 


x  -  |  (1  +  O2 

and  (11) 

y  -  f  n(l  +  O2. 


Point  P  on  the  radial 
where  r  *  -J(x2  +  y^) . 


line  R,  is  at  a  distanee  r  from  the  crack-tip. 
Substituting  from  equations  (11)  yields 


l  „  .  _n2  /, h.2  .  2 

r  -  (1  +  o  y/(j)  +  n 

where 

(1+0 


(12) 


Substituting  equations (11)  in  equation  (7)  yields 


I  (1  +  C)  |  n(l  +  O 

0  J  (1  +  o2 


(13) 


and  its  determinant 
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det | J|  «  jp  (1  +  O3. 

Inverting, 

I  i 

r_l  -1 

11  12 

M  " 

*21  *22 

2 

h(l  +  O 


-4 

h(l  +  €)2 


0  4 

A(1  +  O2 


(14) 


(15) 


Taking  the  derivatives  of  the  displacements  u  and  v  with  respect  to 
€  and  n  gives 


3u 

34 


8 

-  2 
i-1 


3u 

3n 


8 

-  2 

i-1 


3v 

34 


8 

-  2 

i-1 


and 


3v 

3n 


8 

-  2 
i-1 


V 


(16) 


Substituting  for  nodal  numbers  and  performing  the  above  operation  will 
yield 


JV  A  U*  A  U,  A 

■  4~  (-2+3n-n  )  +  u2(l-n)  +  ^(n+n  )  +  jrd-n  ) 

u5  ?  U7  2  u8  2 

-  ^(-n+O  +  u6(l+n)  +  ^(-2-3n-n  )  -  ~  (l-n  )  (17) 

+  (1+  o[^(l-n)-u2(l-n)  +  ^(l-n)  +  ^(l+n)-u6(l+n)+  ^(l+n)l 
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If  -  Cc  +  Ci<l+S), 

and 

-  D„  +  Dl(l+t)  +  Dj  (l+O2 


(20c, d) 


where  A  ,  A, ,  B  ,  B, ,  B„ ,  C  ,  C, ,  D  ,  D, ,  D.  are  constants  for  any 
ololiololz 

given  set  of  nodal  displacements  along  any  line  0  ■  constant. 


Strains  are  equal  to 


"  1“ 

E 

3x 

X 

3v 

£ 

y 

ay 

e 

3u  3v 

xy 

3y  3x 

and 


(21) 


3u  3v 

3u  3v 

3x  3x 

m 

*11  *12 

3£  3C 

3u  3v 

T  T 

3u  3v 

3y  3y 

21  22 
m  • 

3n  3n 

Substituting  equations  (12)  and  (19)  into  equation  (21)  and  performing 


the  operation  will  yield 


1 

1 


(23a-c) 


a 


11 


and 


3v 

3y 


D'  D' 

-2_  +  -l+D' 


(23d) 


where  A^,  A^,  A£,  B^,  Bj,  B^,  C^,  Cj,  C^,  D^,  D^,  D*  are  constants 
for  any  radial  line  (9  “  constant) ,  and  are  independent  of  r. 

If  the  constraints  of  equations  (19a)  and  0.9b)  on  the  nodal 
displacements  at  the  crack-tip  are  imposed,  then  equations  (20a) 
through  (20d)  reduce  to 


3u 

35 


Aq  +  A1(l+5), 


fj-  b1(1+o  +  b2(i+o2. 


(24a-d) 


It  -  co  +  ci(1+e)* 

and 

|s  -  DjU+S)  +  D2(1H)2. 


The  derivatives  will  respect  to  x  and  y  then  reduce  to 


3u 

3x 


3u 

3y 


3v 

3x 


(25a-c) 


12 


and 


3v 

3y 


(25d) 


As  r  -»■  0,  the  terms  In  equations  (23a)  through  (23d)  tend  toward 
(1/r)  while  the  terms  In  equations  (25a)  through  (25d)  tend  toward 
(l//r).  In  both  sets  of  equations,  constant  terms  remain  to  represent 
the  constant-strain  terms  associated  with  thermal  loads.  Therefore, 

If  nodes  1,  7  and  8  at  the  crack-tip  are  restricted  to  having  the 
same  displacements,  the  elastic  inverse  square  root  singularity  will 
occur.  However,  if  nodes,  1,  7  and  8  are  left  free  to  displace 
Independently  of  each  other,  the  perfectly  plastic  Inverse 
singularity  will  occur. 


2.3  BOPACE 

The  finite  element  program  which  was  used  to  model  the  edge- 
cracked  beam  in  tension  was  BOPACE  3D  (The  Boeing  Plastic  Analysis 
Capability  for  Engines)  Version  6.2.  The  program  was  originally 
developed  by  Boeing,  Inc.  for  NASA  for  use  in  the  analysis  of  the 
Space-shuttle  main  engines. 

Many  parameters  which  are  ideal  for  modeling  both  elastic  and 
elastic-plastic  fracture  problems  are  contained  in  the  program.  For 
crack  problems,  a  convenient  option  available  is  the  choice  of 
proportionate  or  serendipity  mapping  [l5] .  It  is  through  the 
serendipity  mapping  that  the  crack-tip  elements  are  formed.  For 
elastic-plastic  analysis,  the  program  employs  an  isotropic- 
kinematic  hardening  theory,  the  Huber-Mises  yield  surface  criterion 

and  the  Prandtl-Reuss  flow  rule  [l9] .  A  choice  of  three  different 
Iterative  schemes  are  also  available  for  the  solution  of  elastic- 
plastic  problems. 


CHAPTER  III 


ANALYSIS  OF  A  SINGLE  EDGE-CRACKED  BEAM 
SUBJECTED  TO  MODE  I  TENSILE  LOADING 

3.1  Method  of  Approach 

The  analysis  of  the  single  edge-cracked  beam  subjected  to  Mode  I 
tensile  loading  was  completed  for  both  the  elastic  and-elastic-perfectly 
plastic  cases.  Two-dimensional  finite  element  models  were  constructed 
for  both  cases.  In  the  elastic  case.  Mode  I  stress  Intensity  factors 
were  numerically  determined  from  crack-tip  nodal  displacements 
[Appendix  Aj .  The  Mode  I  stress  intensity  factors  obtained  from  the 
finite  element  model  nodal  displacements  were  then  compared  to  the 
Mode  I  stress  intensity  factors  obtained  from  a  boundary  collocation 
solution  [Appendix  b].  Both  plane  stress  and  plane  strain  cases 
were  modelled  and  analyzed.  For  the  elastic-perfectly  plastic  cases, 
the  size  and  character  of  the  plastic  zones  were  analyzed. 

3.2  Finite  Element  Meshes 

In  the  elastic  models,  proportionate  elements  (8-node  quads) 
were  used  everywhere  except  for  the  few  elements  immediately 
surrounding  the  crack-tip  where  elastic  crack-tip  elements  were 
employed  (Figures  2  and  3).  Because  of  symmetry,  only  the  upper 
half-plane  of  the  edge-cracked  beam  was  modeled.  Therefore,  all 
nodes  on  the  x-axls  to  the  right  of  and  including  the  crack-tip 
nodes  were  restrained  from  any  displacement  in  the  y-direction.  A 
uniform  tensile  stress  was  applied  to  the  top  row  of  elements.  The 
only  difference  between  finite  element  models  A  and  B  are  the  number 
of  nodes  at  the  crack-tip  and  the  geometry  of  the  crack-tip  elements. 


ENLARGEMENT  OF  CRACK  TIP  ELEMENTS  31-34 


Parameters: 
v  -  .35 
E  -  51.7  x  lo‘ 

“applied  '  601 


Figure  2.  Elastic  Finite  Element  Model  A 


y  f  °APPLIED 


Figure  3.  Elastic  Finite  Element  Model  B. 
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la  the  case  of  the  elastic-perfectly  plastic  models,  proportionate 
elements  (8-node  quads)  were  used  everywhere  except  for  the  few 
elements  immediately  surrounding  the  crack-tip.  Perfectly-plastic 
crack-tip  elements  were  used  to  surround  this  area  (Figures  4  and  5). 
Only  the  upper  half-plane  of  the  edge-cracked  beam  was  modelled. 
Therefore,  all  nodes  on  the  x-axis  to  the  right  of  the  crack-tip 
nodes  were  restrained  from  displacement  in  the  y-direction.  As  in 
the  elastic  case,  a  uniform  tensile  stress  was  applied  to  the  top  row 
of  elements. 

3.3  Elastic  Results  and  Discussion 

The  displacement  results  Immediately  surrounding  the  crack-tip 
of  the  elastic  single  edge-cracked  beam  subjected  to  Mode  I  tensile 
loading  are  presented  in  Figures  6-13.  Finite  element  model  results 
in  terms  of  crack-tip  nodal  displacements  are  compared  with  the 
theoretical  solutions  obtained  from  the  classic  near  field  singularity 
solution  [l]  and  appropriate  stress-intensity  factor.  These 
displacement  results  along  with  additional  results  surrounding  the 
crack-tip  vicinity  are  also  presented  in  tabular  form  in  Appendix  D. 

As  can  be  seen  in  the  figures,  the  best  elastic  results  were 
attained  at  the  quarter-point  nodes  of  finite  element  model  B. 

This  was  the  model  with  9  restrained  nodes  at  the  crack-tip.  The 
difference  between  these  finite  element  model  results  and  the 
theoretical  solution  ranged  from  0-4%  for  both  the  plane  stress  and 
plane  strain  assumptions.  Therefore,  use  of  the  finite  element 
model  displacements  at  the  quarter-point  nodes  would  result  in 
Mode  I  stress  intensity  factors  with  comparable  accuracy. 


Figure  5.  Magnification  of  Crack-tip  Region  in  Elastic-Perfectly 
Plastic  Finite  Element  Model  C. 
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43.0  135.0 


0.0  90.0  180.0 

Angle,  0  (Degrees) 

+  3  Finite  element  model  results  at  quarter-point  nodes 
*  3  Theoretical  solution  results  at  quarter-point  nodes 

Figure  6.  Finite  Element  Model  A  Results  (Plane  Stress  Assumption). 


Displacement  in  x  direction 


Angle,  0  (Degrees) 

Finite  element  model  results  at  quarter-point  nodes 
Theoretical  solution  results  at  quarter-point  nodes 

Finite  Element  Model  A  Results  (Plane  Strain  Assumption) 


Displacement  in  x  direction 
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49.0  199.0 

0.0  90.0  190.0 

Angle.  0  (Degrees) 

+  ■  Finite  element  model  results  at  quarter-point  nodes 
*  ■  Theoretical  solution  results  at  quarter-point  nodes 


Figure  9.  Finite  Element  Model  B  Results  (Plane  Strain  Assumption) 


Angle,  0  (Degrees) 


Finiee  element  model  results  at  quarter-point  nodes 
Theoretical  solution  results  at  quarter-point  nodes 


Finite  Element  Model  A  Results  (Plane  Strain  Assumption) 
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3.4  Elastic-Perfectly  Plastic  Results  and  Discussion 

Results  for  the  elastic-perfectly  plastic  crack-tip  analysis 
are  shown  in  Figures  14-17.  Results  presented  include  examination  of 
the  extent  and  character  of  the  plastic  zone  for  both  plane  stress 
and  plane  strain  conditions  and  predicted  displacements  along  the  free 
flank  of  the  crack.  The  elastic  finite  element  model  crack  flank 
displacements  are  also  presented  to  verify  that  the  model  is  behaving 
correctly  by  comparing  finite  element  displacements  to  those  predicted 
by  the  analytical  elastic  solution. 

In  Figure  14,  plane  stress  was  assumed  for  the  two-dimensional 
finite  element  model  and  the  characteristic  plane  stress  plastic 
zone  shape  was  attained  [2l] .  Immediately  outside  the  plastic  zone, 
the  y-displacements  were  examined  along  the  crack  flank  to  determine 
whether  they  still  followed  a  1//F  relationship.  As  is  evident,  the 
y-displacements  from  the  elastic-perfectly  plastic  finite  element 
model  were  greater  in  magnitude  than  both  the  predicted  elastic  values 
and  the  elastic  finite  element  values.  Not  only  were  the  y-displacements 
greater  in  magnitude,  but  they  did  not  follow  the  \/\fc  relationship 
characteristic  of  the  elastic  displacements.  This  can  be  explained  by 
the  increased  compliance  due  to  the  plastic  zone.  Another  approach 
was  therefore  attempted.  In  this  approach,  the  concept  of  an 
equivalent  crack  length  was  introduced.  In  order  to  predict  the 
elastic-perfectly  plastic  results,  the  same  linear  elastic  Mode  I 
displacement  equations  are  used  with  the  plastic  zone  being  treated 
as  nothing  more  than  an  extension  to  the  crack  length,  ry.  The  origin 
of  the  axis  system,  however,  is  still  located  at  the  original  crack- 
tip.  A  similar  approach  was  adopted  in  the  past  by  F.  A.  McClintock 


Plane  Stress  Plastic  Zone  Observed  in  Elastic-Perfectly 
Plastic  Finite  Element  Model  C. 


'-T-T-T^T-T-T^T-T-VT'W-V-VT'IV-' 


Figure  15.  Plane  Strain  Plastic  Zone  Observed  in  Elastic-Perfectly 
Plastic  Finite  Element  Model  C. 


y-displacements,  v  x  10  (inches) 

E  *  Elastic  finite  element  model  results  along  crack  flank 
P  =  Elastic-perfectly  plastic  finite  element  model  results 
along  crack  flank 

+  =  Predicted  elastic  y-displacement  along  crack  flank 
*  *  Predicted  y-displacements  along  crack  flank  using  the 
equivalent  crack  length 

Figure  16.  Elastic/plastic  Results  (Plane  Stress  Assumption). 


y-displacement,  v  x  10  (inches) 

E  =  Elastic  finite  element  model  results  along  crack  flank 
P  =  Elastic-perfectly  plastic  finite  element  model  results 
along  crack  flank 

+  =  Predicted  elastic  y-displacements  along  crack  flank 
*  *  Predicted  y-displacements  along  crack  flank  using  the 
equivalent  crack  length 

Figure  17.  Elastic/Plastic  Results  (Plane  Strain  Assumption). 
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and  G.  R.  Irwin  In  an  analysis  for  the  Mode  III  case  [22] .  Figure  16 
shows  the  results  from  this  approach.  Again,  the  predicted  y-displace- 
ments  attained  by  using  the  equivalent  crack  length  did  not  fit  the 
elastlc-perfectly  plastic  finite  element  model  results  but  did 
approach  them  asymptotically  with  increasing  distance  from  the  crack-tip. 

In  Figure  15,  plane  strain  was  assumed  for  the  two-dimensional 
finite  element  model  and  the  characteristic  plane  strain  plastic 
zone  was  attained  [2l] .  Again  as  in  the  plane  stress  model,  the 
y-displacements  along  the  crack  flank  outside  the  plastic  zone  were 
examined  to  see  if  they  followed  a  1//F  relationship.  However,  the 
y-displacements  from  the  elastic-perfectly  plastic  finite  element 
model  were  greater  in  magnitude  than  the  predicted  elastic  values 
and  did  not  show  the  l//r  relationship  (Figure  17) .  As  in  the  plane 
stress  case,  the  equivalent  crack  length  approach  was  again  used  to 
predict  the  plane  strain  elastic-perfectly  plastic  results.  These 
results  are  presented  in  Figure  17.  Again  the  predicted  y-displace- 
ments  attained  by  using  the  equivalent  crack  length  did  not  fit  the 
elastic-perfectly  plastic  finite  element  model  results  but  did 
approach  them  asymptotically  with  increasing  distance  from  the  crack- 


tip. 


CHAPTER  IV 


SUMMARY  AND  CONCLUSIONS 

By  constructing  various  two-dimensional  finite  element  models 
for  a  single  edge-cracked  beam  subjected  to  Mode  I  tensile  loading, 
a  thorough  analysis  of  the  stresses,  strains  and  displacements 
surrounding  the  crack-tip  was  performed.  Particular  emphasis  was 
placed  on  1)  determining  the  Mode  I  stress  intensity  factor  in  the 
elastic  models  via  a  displacement  method  and  2)  examining  the  size 
and  character  of  the  plastic  zone  in  the  elastic-perfectly  plastic 
models.  In  all  of  the  finite  element  models,  either  special  elastic 
or  plastic  crack-tip  elements  were  used  to  create  the  crack-tip 
singularity. 

The  elastic  finite  element  model  results  do  indeed  represent 
accurate  crack-tip  singularity  effects  as  was  made  evident  by  a 
comparison  with  available  analytical  solutions.  Not  only  do  the 
elastic  models  represent  the  singularity  effect  accurately,  but  with 
the  aid  of  the  special  elastic  crack-tip  elements,  they  are  able  to 
do  it  with  a  minimal  amount  of  elements.  This  leads  to  a  significant 
savings  in  both  time  and  costs 

The  elastic-perfectly  plastic  finite  element  model  results  were 
also  of  great  interest.  The  plastic  zones  obtained  from  the  finite 
element  models  did  indeed  agree  with  the  characteristic  plane  stress  and 
plane  strain  plastic  zones.  However,  the  equivalent  crack  length 
approach  did  not  yield  accurate  results  for  the  y-displacements  along 
the  crack  flank  for  either  the  plane  strain  or  plane  stress  case. 


CHAPTER  V 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

I 

A  great  deal  of  additional  research  can  and  should  be  conducted 
in  applying  the  finite  element  method  to  fracture  mechanics  problems. 
Uith  the  improved  finite  element  programs  available,  several  promising 
areas  can  be  investigated. 

One  area  which  can  be  pursued  further  is  the  investigation  of  a 
three-dimensional  elastic  finite  element  model  for  a  single  edge- 
cracked  beam  subjected  to  Mode  I  tensile  loading.  With  a  three- 
dimensional  model,  one  could  determine  the  effect  of  model  thickness 
on  stress  state  (plane  stress  or  plane  strain)  at  the  crack-tip. 

Finite  element  model  results  could  then  be  compared  to  empirical 
and  experimental  results.  One  experimental  technique  which  has 
been  used  to  investigate  this  area  is  scattered-light  photo- 
elasticity  [20]  . 

Another  area  which  should  be  pursued  is  the  investigation  of  a 
three-dimensional  elastic-perfectly  plastic  finite  element  model  for  a 
single  edge-cracked  beam  subjected  to  Mode  I  tensile  loading.  With 
a  three-dimensional  elastic-perfectly  plastic  finite  element  model 
a  three-dimensional  characterization  of  the  plastic  zone  surrounding 
the  crack-tip  line  could  be  attained. 

Finally  a  thorough  Investigation  of  the  J-integral  as  it  is 
applied  in  the  area  of  elastlc/plastlc  fracture  mechanics  should  be 
considered . 

In  any  investigation,  a  careful,  well  drawn  out  in  advance, 
stepwise  procedure  from  the  simple  to  the  more  complex  case  should 
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be  undertaken.  It  la  this  type  of  procedure  which  allows  one  to 
gain  valuable  experience  with  the  actual  "inner  workings"  of  every 
aspect  of  the  computer  model  and  will  assure  one  of  accurate 
results.  Therefore,  in  the  investigation  of  three-dimensional 
finite  element  models,  simple  models  with  coarse-grids  should  be 
constructed  first.  These  models  should  then  be  thoroughly 
analyzed  to  verify  that  everything  is  operating  correctly.  If  any 
problems  do  occur,  they  are  easily  seen  and  can  be  corrected. 

Once  everything  has  been  checked  and  is  operating  correctly  in  the 
simple,  coarse-grid  models,  larger  and/or  finer-grid  models  may  be 
generated. 
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APPENDIX  A 

ELASTIC  CRACK-TIP  STRESS  AND  DISPLACEMENT  FIELDS 


Refering  to  Figure  Al,  the  well-known  Mode  I  stress  and 
displacement  fields  surrounding  a  crack-tip  [l,  23]  are  given  by 
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where  F  -  for  the  plane  stress  assumption, 

and  F  *  3-4v  for  the  plane  strain  assumption. 

APPENDIX  B 


SINGLE  EDGE  NOTCH  TEST  SPECIMEN 


From  reference  14  the  Mode  I  stress  intensity  factor  for  the 
single  edge-cracked  plate  (Figure  Bl)  is  given  as: 


Kj  *  o  /iia  F(a/b) , 


Empirical  formulas  for  F(a/b)  are  given  as: 

1.  F(a/b)  -  1.12  -  .231  (a/b)  +  10.55  (a/b)2  -  21.72  (a/b)3 


+  30.39  (a/b) 

Gross  1964,  Broun  1966;  Boundary  Collocation  Method 

Accuracy:  .5Z  for  a/b  £  .6 

2.  F(a/b)  -  .265  (1-a/b)4  +  265  <a/b) 

(l-a/b)3/2 

Tada  1973 

Accuracy:  1Z  for  a/b  <  .2;  .5Z  for  a/b  £  .2 
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Accuracy:  .5%  for  any  a/b 
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APPENDIX  C 

STRAIN  ENERGY  RELATIONS 

From  strain  energy  relations  the  strain  energy  release  rate  is 
given  as 


v 

G*  - —  for  plane  stress, 

X  £ 

(C-l) 

and 
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G  “  — ■  v—  K_^  for  plane  strain. 

X  E  X 

(C-2) 

Rearranging  equations  (C-l)  and  (C-2)  yields 

Kj.  ■  v/eGj.  for  plane  stress 
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and 
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